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Abstract 

We propose a test for non-standard neutrino-neutrino interactions by using ul- 
trahigh energy AGN neutrinos. Such interactions would influence the AGN 
neutrino flux due to collisions with cosmic background neutrinos. For typical 
AGN neutrinos we obtain an upper limit for the coupling constant g < 6.4- 
if the mediator is light and g/(Mx/MeV) < 0.013 if the mediator is heavy. We 
compare our results with constraints from other phenomena previously consid- 
ered. 
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1. Introduction. Active galactic nuclei (AGN) are found to be the origin of the 
highest energy photons detected in gamma ray observatories. The production of these 
photons is believed to be accompanied by a production of other high-energy cosmic 
rays, like protons, electrons and neutrinos ffl. The weakly interacting neutrinos can 
carry us first-hand information of the processes taking place in the source and of the 
source itself, and they can reach us from more distant AGN. On the other hand, 
the AGN neutrinos can be useful in studying the properties of neutrinos themselves. 
The high-energies not achievable in laboratory experiments or in other astrophysical 
environments, combined with long distances between the source and the detection 
point, may reveal new important information of mixing, masses and possible new 
interaction forms of neutrinos. 

In this paper we will study the possibility of testing non-standard neutrino-neutrino 
interactions with high-energy neutrinos from AGN. If such interactions exist, the AGN 
neutrinos will undergo extra collision processes with neutrinos in cosmic background 
in their way from the source to detectors. This will affect their measured flux. We 
will derive an upper limit for the effective strength of these interactions that could be 
reached by measuring the flux in the new neutrino telescopes such as AMANDA 0, 
Nestor ||, Baikal and ANTARES f|. The energies of the extragalactic neutrinos 
we are interested in are typically from some hundreds of GeV up to a few EeV:s. 
The typical length scale is several hundreds of megaparsecs which is the distance to 
the neighbouring active galaxies. The new neutrino telescopes can detect very high 
energy neutrinos with quite a good angular resolution, in the future even allowing 
one to connect observed neutrinos with given AGN sources despite the atmospheric 
neutrino background ||. 

In the standard model neutrinos interact with other neutrinos only via Z-boson 
exchange. Neutrino- neutrino interactions are, however, hard to test, and it is not 
excluded that there are substantial non-standard "secret" forces between neutrinos. 
We will consider the possibility of testing such secret forces using AGN neutrinos. To 
be specific, we will assume that these interactions are of the form {i = e, fi, r) 

gva^iX 1 ", (1) 

where A is a spin-one boson and g is the coupling constant of the interaction. We will 
here assume that the coupling is similar between different neutrinos and that this new 
vectorial interaction does not couple X boson with any other particles than neutrinos. 
The effective Hamiltonian of this interaction in the low energy approximation is 

H = F v (u la u)(u 1 a u), (2) 

where Fy = g 2 /Mjr is the effective coupling constant of the interaction. In the order 
of magnitude level our results are obviously valid for a wider class of models than just 
for the one defined in Eq. ([]]) and Eq. (Q). 



Laboratory constraints for the non-standard interactions of the form given above 
have been discussed in - [10]. The most stringent one comes from measuring the 
invisible width of the Z-boson at LEP [[HJ: 



F v < 4 x 10 2 G F . (3) 

Here the number of neutrino species have been used to set a limit on the non-standard 
four neutrino coupling resulting in the cascade decay Z — > vv — > vvvv. 

Constraints from cosmology and astrophysics have been considered in [[Ll|] - [ |i~3|| . 
One of these limits [[H| is even stricter than the laboratory limit (H), but it is more 
model dependent, as will be discussed below. These constraints come from SN1987A JIT 



12] and nucleosynthesis considerations [13 



It is interesting to study whether one could improve the present bounds by mea- 
suring the flux of AGN neutrinos. We will show that indeed for certain values of the 
X boson mass more stringent limits will be achievable. 

2. UHE neutrino sources. In the litterature two types of mechanisms for neutrino 
production in AGN have been discussed, spherical accretion mechanism and blazars 
producing jets (for a review, see [0]). 

In the spherical accretion mechanism it is believed that in the center of an active 
galaxy a supermassive black hole creates an accretion disk of materia falling into the 
black hole. In this disk a shock is formed that accelerates protons [14], [15] to very 
high energies. Neutrinos are produced in pion photoproduction processes and in pp- 
collisions, the most important reaction chain being pry — > mr + — > — > e + i/ e u fJl [ [U| . 



Different models have been presented |17], [18], [19], ^1 f° r estimating the fluxes of 



nucleons and neutrinos starting from the observed X ray flux. 

It is argued that in blazars protons are accelerated in jets [^1 - [22] by the first 



order shock acceleration mechanism, and neutrinos are then created in the decays 
of photoproduced pions. The flux of neutrinos can be calculated from the observed 
photon spectrum |23|| . 

In both cases of AGN tau neutrinos are produced only negligible amounts, but 
they might appear due to mixing. In case the secret interactions exist, the possible 
flux of tau neutrinos will be reduced by the same amount as that of u e and v^. 

The detected flux of neutrinos from a single source is estimated to be from a few 
to even some hundreds of events per year in a square kilometer neutrino detector |J. 
However, the statistics will improve with time and when the plans to build several 
new kilometer scale detectors are to come reality. Also, the higher is the neutrino 
energy, the better is its detection rate in a neutrino telescope compared with that of 
the background atmospheric neutrinos. Moreover, one can increase the statistics by 
integrating the flux of all AGN and blazars of the whole known Universe. The inte- 



grated neutrino flux is expected to be well above the background flux of atmospheric 
neutrinos |]. 

To obtain a conservative limit on the strength of the secret interactions we ap- 
proximate the distance of AGN and blazar sources to be 500 Mpc, which is less than 
the actual distance of the majority of them. More stringent limits will be obtained 
by statistically indentifying the individual sources of the observed neutrinos and using 
the corresponding actual distance. 

3. An order of magnitude estimate. Constraints on the secret neutrino-neutrino in- 
teractions will be obtained from the possible depletion of the AGN neutrino flux. To 
obtain an order of magnitude estimate we will compare the mean free path A of AGN 
neutrinos associated with these interactions and the distance D to the source with each 
other. An upper limit for the strength of the force is set by the condition \~ 1 D < 1. 

Assuming that the background neutrinos are relativistic Dirac particles and that 
the interaction is of the vector type (|1|), one can estimate the cross section of an 
ultrarelativistic neutrino colliding with a cosmic background neutrino as 



where s cx 2E]T, E\ is the energy of the AGN neutrino and T = 1.9 K = 1.64- 10 -4 eV 
is the temperature of the background neutrino. Here we have not taken into account 
the thermal distribution of the neutrino background. If ^/s ^> Mx, then the mass of 
the mediator can be neglected and the cross section of Eq. can be approximated 
as a = g 4 /s. On the other hand, if ^/s <C Mx, then the cross section can be written 
as o = g A s/M x . 

The mean free path of the high energy neutrino is A = l/(n v a), where n v is the 
number density of the background neutrinos, which for one neutrino family is n Vi = 
n 9i = 55 cm -3 . For representative values of incoming neutrino energy of E% = 10 15 eV 
and distance to the source of D = 500 Mpc = 1.5 • 10 27 cm, the condition \~ X D < 1 
leads, assuming the neutrino flux to be according to the theoretical models, for the 
coupling constant an upper bound of 



g < {s/{n v \)) 1/A ~ 0.01 



(5) 



in the case \fs ^> Mx and 



#/(M x /MeV) < {n u Xs)- 1/4 ~ 0.02 



(6) 



in the case s/s M x . 

In the next section we will define this crude estimate by adding up the cross sec- 
tions of all relevant processes and also taking into account the thermal distribution 



of the background neutrinos. A similar analyses for supernova SN1987A neutrinos is 
presented by Kolb and Turner in [12 . 



/^/(p 2 )|Av|„W, (7) 



4- Thermal analysis. We can write the Boltzmann equation of high energy neutri- 
nos Vi with phase space distribution /' and background neutrinos u 2 , assumed to be 
relativistic, with Fermi-Dirac distribution / in the form 

f dt J (2tt) 

where s = (p± +P2) 2 and Av = v 2 — vi. Here we have omitted the usual Pauli blocking 
factors as there is no Fermi degeneracy in the cosmic background. The cross sections 
a(s) for each different process are given in the table [l] below, where we have taken as 
an example high energy electron neutrino scatterings with background neutrinos or 
antineutrinos and annihilations producing neutrino-antineutrino pairs. The numerical 
values of the cross sections on the last three rows are multiplied with two since there 
are two possible final states. Muon and tau neutrinos and antineutrinos undergo the 
corresponding reactions. 

In terms of y = i|vi|, where t is the flight time from the source to the detector and 
y the distance, the Boltzmann equation reads 

The right-hand side of Eq. (H) is the inverse of the mean free path, A -1 . We have here 



approximated that s ~ 2E 1 E 2 (1 — z) and |Av| ~ y 2(1 + z), where z = cos# and 6 is 
the angle between the initial state momenta. 

The total cross section of the reactions an AGN neutrino undergoes is of the form 
a = a/s and a = as/M x in the case \fs 3> Mx and y/s <C Mj, where the dimension- 
less constant a will be calculated numerically below in table [l|. The mean free path of 
the AGN neutrinos is given by the expression 



, a J °°dE 2 E 2 f(E 2 ) _a ^ 

E x Jo 00 dE 2 E 2 f(E 2 ) 712 E\T 12 C(3) ' 



^ = ^ roo:jJ 2£ )J> 2 = ^77777^^ (9) 
'-1 - 

in the case \fs ^> M x and 

t lQ aE 1 J™dE 2 Elf(E 2 ) !QaE 1 T 7n 4 
' ~ h M x ^dE 2 Elf{E 2 ) n2 5M1 180C(3) n2 * 



in the case \fs Mx, where n 2 is the number density of the cosmic background 
neutrinos. 

As can be seen from the table 0, some of the cross sections are divergent when 
t = — s or t = 0. In the case the cross section has a divergency at t = —s, the 
integration is started from t = — (1 — e)s where e is set to be 0.1. If there is a 
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2 x 0.053 




s 2 ,(s + t) 2 
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2 x 0.569 


(s+t) 2 +s 2 


2 x 0.053 




t 2 ,(s + t) 2 

7* + s 2 


2 x 0.027 


(s+t) 2 +t 2 


2 x 0.027 


constant a 




4.32# 4 




0.558# 4 



Table 1: Processes and cross sections. The differential and integrated numerical values 
of the cross sections of the processes are listed in the second and third column in the 
case ^/s ^> Mx and in the last two columns in the case \fs <C Mx- 

divergency at t — 0, the integration is performed up to t = —es. This procedure will 
not substantially affect the flux of high energy neutrinos: when t — * 0, the collision 
is elastic and no reaction occurs, and when t —>■ —s, the final state neutrino has the 
same momentum as the initial state neutrino. 

Requiring now \~ l D < 1, where D is the distance to the source we end up with 
constraints well accordance with our crude estimates given in Eq. (|^) and Eq. (Q): 

g < 6.4. io-3 ( (f^tyy an 

y V £>/(500 M Pc) 

in the case y/s ^> Mx and 

9 /(M x /UeV) < 0.013 ( (£i/PeV)( ^ /500 Mpc) )' /4 (12) 

in the case yfs -C Mx- If the secret interactions do exist, then depending on the mass 
of the mediator either the low energy (^/i ^> Mx) or the high energy (a/s -C Mx) 
part of the AGN neutrino spectrum will be influenced. 

5. Comparison with other constraints. Let us make a comparison of our constraints 
with the bounds previously obtained from other phenomena. In references [0], 0, |Q 
upper bounds to secret neutrino-neutrino interactions have been derived by using light 
meson decays. In the LEP constraint ([|), based on the measurement of the invisible 
width of the Z boson, the mediator of the secret interaction is assumed to be a vector 
particle much heavier than the Z boson, i.e. My 3> Mz- In order to compare this 
bound with the one we obtained, we must extrapolate the LEP constraint to lower 
mediator masses. The coupling constant is of the form Fy = g 2 /M v < 400Gp, and if 
we set My = 100 GeV, we will get g 2 < 4.0 x 10 6 Gf- A comparison with our result (p~2j) 
with representative values E\ — 1 Pev and D = 500 Mpc, giving g/Mx < 13/GeV, 
shows that our constraint is more stringent whenever Mx < 0.5 GeV. 



From supernova SN1987A data two kinds of limits were obtained to the neutrino- 



neutrino couplings, considered in [I ![ and 1 12| . In [I l[ the constraint was derived 



from the duration of the neutrino pulse and the magnitude of neutrino flux. A limit 
to the neutrino-neutrino cross section was placed to cr pp < 1CT 35 cm 2 , and g 2 < 
5.1 • 10~ 6 (T I/ /10 MeV)rj It was assumed that the possible secret interaction was due 
to the exchange of a scalar particle with mass smaller than the neutrino temperature 
T u , which is less than 10 MeV. Our result is of the same order of magnitude when this 
limit is interpreted in terms of vector particle exchange. 

In [12] the supernova neutrino interactions with cosmic neutrino background was 
used to derive limits for the couplings. The energy of the supernova neutrinos is 
around 10 MeV and the distance 55 kpc is much shorter than that of AGN. The 
limits obtained are comparable with those we obtained, g < 5.6 • 10 -4 if My <C 60 eV 
and g/(M v /MeV) < 12 if My > 60 eV (in this case y/s ~ 60 eV). Comparing these 
results with our equations ([□]) and (|i~2D shows that the supernova constraints are more 
stringent than the AGN constraints in the case the mass of the mediator is below 1 
keV (assuming again E\ — 1 PeV and D = 500 Mpc). With larger mass values the 
AGN limits always surpass the supernova ones. For example with Mx = 10 keV we 
obtain g < 6.4 • 10~ 3 whereas the supernova limit is g < 0.12. 

The nucleosynthesis limit was obtained [|13|] by requiring that right-handed Dirac 
neutrinos have to decouple at temperatures above QCD phase transition. This leads 
to an upper limit for the coupling strength which in the case of massive mediator is 
obtained to be Fy < 3 x 10~ 3 GV and in the case of massless mediator g < 2 x 10 -5 . 
These limits are actually more stringent than any other limit discussed above, includ- 
ing ours, but they are, however, valid only for a very limited class of models. 



6. Summary. We have considered a testing of non-standard neutrino-neutrino inter- 
actions by using ultrahigh energy AGN neutrinos. If such interactions exist, AGN 
neutrino flux should be influenced by extra collisions with cosmic background neu- 
trinos. We have derived the upper limit for the coupling strength one can probe by 
measuring the AGN neutrino flux. Depending on the mass of the mediator of these 
secret interactions the constraints we obtained, given in equations (|TT|) and (|T2|), are 
more stringent than previously presented astrophysical and laboratory limits. 
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1 According to reference [^4| this result is disputable and the upper limit could actually be some- 
what larger. 
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